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Aan  schokreductiemaatregelen,  met  name  het  verend  opstellen  van  apparatuur 
wordt  vanzelf sprekend,  al  sinds  lange  tijd  veel  aandacht  en  onderzoek 
besteed. 

Met  name  de  laatste  jaren  is  het  belang  van  schokreductie  nog  verder 
toegenomen,  onder  invloed  van  krimpende  defensiebudgetten  en  de  daarmee 
gepaard  gaande  wens  en  verwachting  om  bij  nieuwbouw  tot  kostenbesparing  te 
geraken  door  ondermeer  specifieke  schokeisen  aan  apparatuur  naar 
toeleveranciers  waar  mogelijk  te  laten  vervallen. 

Onder  deze  omstandigheden  is  het  van  belang  dat  verende  opstellingen  een  zo 
groot  mogelijke  schokreductie  bieden  zodat  het  stellen  van  schokeisen  aan 
apparatuur  dan  inderdaad  zonder  al  te  groot  risico  achterwege  kan  blijven. 
Omdat,  meer  nog  dan  vroeger,  zeer  grote  slap  verend  opgestelde  vloeren 
zullen  worden  gaan  toegepast,  verdienen  deze  extra  aandacht. 

Eerder  uitgevoerde  schokresponsieberekeningen,  [1],  [2],  [3],  hadden 
betrekking  op  representatieve  dek  -  vloer  combinaties  zoals  mogelijk 
aanwezig  in  de  commandocentrale  van  het  LCF. 

Bij  de  laatste  van  deze  berekeningen  [3]  was  sprake  van  een  vloer  opgebouwd 
uit  door  de  KSG  voorgestelde  stalen  UNP  80*50*5  profielen  en  opgesteld  op 
kabelgewonden  veren  van  Sebert .  Van  die  dek-vloer  combinatie  werden 
vervolgens  eigentrillingsyormen  en  frequenties  bepaald,  maar  ook  van  het 
dek  en  de  vloer  afzonderlijk  en  tevens  voor  enkele  gemodif iceerde 
constructies  [4] . 

Daarmee  werd  het  inzicht  in  de  interaktie  tussen  dek  en  vloer  verder 
vergroot . 

Uit  die  yorige  berekeningen  [3] ,  [4]  kwam  onder  meer  naar  voren  dat  de 
doorgeleidde  versnelling  naar  de  vloer  al  oscillerend  aanzienlijk  toenam  en 
groter  werd  dan  verwacht. 

Daarom  werd  nu  door  middel  van  een  nieuwe  schokresponsieberekening  verder 
aandacht  geschonken  aan  het  effect  van  het  vrijwel  samenvallen  van  de 
eigenf requenties  van  dek-  en  vloertrilvormen. 

Tot  dusverre  was  in  de  rekenmodellen  alleen  de  demping  verdisconteerd  van 
de  kabelgewonden  veren.  De  onbekende  demping  in  de  constructie  was  niet 
meegenomen.  In  een  tweede  nieuwe  responsieberekening  is  nu  demping  in  het 
dek  ingebouwd. 
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SAMENVATTING 

Aan  schokreductiemaatregelen,  met  name  het  verend  opstellen  van  apparatuur 
wordt  vanzelfsprekend,  al  sinds  lange  tijd  veel  aandacht  en  onderzoek 
besteed. 

Met  name  de  laatste  jaren  is  het  belang  van  schokreductie  nog  vender 
toegenomen,  onder  invloed  van  krimpende  defensiebudgetten  en  de  daarmee 
gepaard  gaande  wens  en  verwachting  om  bij  nieuwbouw  tot  kostenbesparing  te 
geraken  door  ondermeer  specifieke  schokeisen  aan  apparatuur  naar 
toeleveranciers  waar  mogelijk  te  laten  vervallen. 

Onder  deze  omstandigheden  is  het  van  belang  dat  verende  opstellingen  een  zo 
groot  mogelijke  schokreductie  bieden  zodat  het  stellen  van  schokeisen  aan 
apparatuur  dan  inderdaad  zonder  al  te  groot  risico  achterwege  kan  blijven. 
Omdat,  meer  nog  dan  vroeger,  zeer  grote  slap  verend  opgestelde  vloeren 
zullen  worden  gaan  toegepast,  verdienen  deze  extra  aandacht. 

Eerder  uitgevoerde  schokresponsieberekeningen,  [1],  [2],  [3],  hadden 

betrekking  op  representatieve  dek  -  vloer  combinaties  zoals  mogelijk 
aanwezig  in  de  commandocentrale  van  het  LCF. 

Bij  de  laatste  van  deze  berekeningen  [3]  was  sprake  van  een  vloer  opgebouwd 
uit  door  de  KSG  voorgestelde  stalen  UNP  80*50*5  profielen  en  opgesteld  op 
kabelgewonden  veren  van  Sebert .  Van  die  dek -vloer  combinatie  werden 
vervolgens  eigentrillingsvormen  en  frequenties  bepaald,  maar  ook  van  het 
dek  en  de  vloer  afzonderlijk  en  tevens  voor  enkele  gemodif iceerde 
constructies  [4] . 

Daarmee  werd  het  inzicht  in  de  interaktie  tussen  dek  en  vloer  verder 
vergroot . 

Uit  die  vorige  berekeningen  [3] ,  [4]  kwam  onder  meer  naar  voren  dat  de 

doorgeleidde  versnelling  naar  de  vloer  al  oscillerend  aanzienlijk  toenam  en 
groter  werd  dan  verwacht. 

Daarom  werd  nu  door  middel  van  een  nieuwe  schokresponsieberekening  verder 
aandacht  geschonken  aan  het  effect  van  het  vrijwel  samenvallen  van  de 
eigenfrequenties  van  dek-  en  vloertrilvormen. 

Tot  dusverre  was  in  de  rekenmodellen  alleen  de  damping  verdisconteerd  van 
de  kabelgewonden  veren.  De  onbekende  damping  in  de  constructie  was  niet 
meegenomen.  In  een  tweede  nieuwe  responsieberekening  is  nu  damping  in  het 
dek  ingebouwd. 
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1 .  INLEIDING 


In  1995  ward  een  meer jarenproject  "Schokreductiemaatregelen"  opgestart. 

In  het  kader  van  dat  project,  dat  niet  exclusief  gericht  is  op  het  LCF, 
wordt  onder  meer  aandacht  bestead  aan  de  schokresponsie  van  grote  zwevende 
vloeren.  Om  daarbij  aan  te  sluiten  bij  reele  afmetingen  en  constructies 
ward  op  verzoek  van  de  Af deling  Scheepsbouw  gerekend  aan  een 
representatieve  dwarsscheepse  doorsnede  van  een  dek  en  de  daarop  slap 
verend  opgestelde  vloer  t.p.v.  de  commandocentrale  van  het  LCF. 

In  het  kader  van  dat  project  zijn  daarover  tot  nu  toe  dx'ie  interim 
rapporten  verschenen:  [1] ,  [2]  en  [4],  terwijl  in  het  kader  van  een  ander 
project,  dat  wel  direct  betrekking  heeft  op  het  LCF,  nog  een  vierde  rapport 
verschenen  is  [3] . 

In  het  kader  van  dat  meerjarenproject  "Schokreductiemaatregelen"  zal 
blijkens  het  overeengekomen  projectplan  [5]  onder  meer  nog  de  invloed  van 
whipping  onderzocht  worden. 

Kort  samengevat  is  in  de  vier  eerder  uitgebrachte  rapporten  sprake  van  het 
volgende : 

Het  betreft  steeds  een  2D  dwarsscheeps  rekenmodel  met  een  dekbreedte  van 
16.8  m  en  een  vloerbreedte  van  13.2  m. 

De  mate  van  detaillering  in  het  eindige  elementen  model  is  gering,  met  in 
totaal  slechts  een  52-tal  balkelementen  en  een  21-tal  veerelementen. 

Het  eerst  verschenen  rapport  [1]  betreft  de  schokresponsie  van  een  dek- 
vloer  combinatie  bij  toepassing  van  kabelgewonden  veren  van  het  type 
Polycal  Y1642-2  van  Socitec,  waarbij  de  vloer  is  opgebouwd  uit  aluminium 
pallet  delen.  Qua  aanstoting  is  in  de  eerste  plaats  de  gestandaardiseerde 
schokbeweging  pulsvorm  C(F)  toegepast  (bij  som  3)  en  vervolgens  de  kickoff 
puls  (bij  som  4) .  Geconcludeerd  werd  dat  deze  laatste  het  meest  realistisch 
is . 

Het  tweede  rapport  [2]  betreft  de  zelfde  dek-vloer  combinatie,  maar  nu  met 
toepassing  van  rubber  veren  van  het  type  N2C  kl.4  van  Loggers. 

Deze  veren  worden  ook  wel  toegepast  onder  dit  soort  vloeren  in  Duitse 
f^G^stten.  Deze  laatste  berekeningen  (som  6)  lieten  grotere  doorgeleide 
versnellingen  naar  de  vloer  zien  (  200  m/s^'  vergeleken  met  60  m/s^  bij  de 
kabelgewonden  veren)  en  ook  grotere  buigspanningen  in  de  vloerbalk. 

Het  derde  rapport  [3]  heeft  betrekking  op  een  geheel  andere  vloer.  Onlangs 
heeft  de  KSG  namenlijk  een  tekening  uitgebracht,  getiteld  "Voorstel 
Schokvloeren" ,  gedateerd  31-5-96,  betrekking  hebbend  op  de  computerruimte 
Sewaco  achter.  Dat  betekent  vervanging  van  de  aluminium  vloer  met  aan 
elkaar  geboute  pallets  door  een  enkel  gelast  stalen  frame,  opgebouwd  uit 
koudgewalste  UNP  80*50*5  profielen.  Deze  zelfde  profielen  zijn  nu  aanwezig 
verondersteld  in  de  vloer  t.p.v.  de  commmandocentrale .  De  vloer  is  wederom 
opgesteld  op  kabel  gewonden  veren,  nu  van  het  type  K19  100-02  van  Sebert . 
Aanstoting  wederom  met  de  kickoff  puls. 

De  uitgevoerde  schokresponsieberekening  wordt  in  [3]  aangeduid  met  som  7. 
Daarbij  ligt  een  vergelijking  voor  de  hand  met  de  resultaten  van  som  4 
omdat  bij  beide  berekeningen  sprake  is  van  aanstoting  met  de  kickoff  puls 
en  ook  van,  qua  veerkarakteristiek,  geli jkwaardige  schokveren. 

Opvallend  is  dan  dat  de  slappere  stalen  vloer  bij  het  meest  linkse  kabinet 
1  een  flinke  opslingering  vertoont  waarbij  de  versnelling  lokaal  aangroeit 
tot  een  waarde  van  100  m/s2,  hetgeen  hoog  is  in  vergelijking  met  de  60  m/s^ 
bij  som  4.  Verondersteld  werd  ([3], biz  13)  dat  er  "blijkbaar  sprake  is  van 
een  eigentrillingsvorm  van  de  vloer  die  dicht  ligt  bij  de  excitatie 
frequentie  van  33  Hz  van  het  dekveld. 

Dit  werd  later  bevestigd  bij  de  berekening  van  de  trilvormen  (som  8) .  Er 
blijkt  eentrilvorm  met  een  frequentie  van  32.7  Hz  te  bestaan  (zie  [4],  fig 
12)  waarbij  vrijwel  alleen  het  dek  en  met  name  daarvan  de  zijvelden 
trillen.  Vlak  daarbij  in  de  buurt  ligt  een  trilvorm  ([4],  fig  11)  met  een 
eigenf requentie  van  31.7  Hz,  waarbij  met  name  de  vloer  trilt  en  het  dek 
vrijwel  stilstaat.  Die  trilvorm  heeft  zijn  grootste  amplituden  bij  het 
linker  kabinet  1 . 

Een  duidelijk  voorbeeld  van  zo'n  lokale  opslingering  van  de  vloer  betreft 
b.v.  knooppunt  5  waarvoor  de  versnelling  weergegeven  werd  gedurende  500  ms 
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(soiti  7,  [3],  fig  A. 7.6).  Deze  figuur  wekt  de  indruk  dat  bij  ongeveer  500  ms 

de  maximal©  amplitude  wel  ongeveer  bereikt  is. Dat  zou  wijzen  op  een  zweving 
met  een  frequentie  van  ongeveer  1  Hz.  Dat  is  inderdaad  het  verschil  in 
eigenf requentie  van  de  twee  genoemde  trilvormen. 

Bij  de  zelfde  som  slingert  b.v.  knooppunt  27  (midden  onder  kabinet  2)  in 
het  geheel  niet  op  ([3],  fig  A. 7. 9).  Dat  is  niet  verwonderli jk  omdat  een 
van  detwee  trilvormen  daar  juist  een  knoop  vertoont  {[4]  , fig  11). 

Die  trilvorm,  met  de  frequentie  van  31.7  Hz,  doet  vermoeden  dat  de  vloer 
ook  onder  het  rechter  kabinet  3  een  aanzienlijke  aangroei  van  de  amplitude 
te  zien  zal  geven  waarbij  de  knooppunten  37  en  49  dan  in  tegenfase  moeten 
bewegen.  Dat  is  inderdaad  het  geval  blijkens  de  in  het  voorliggende  rappoi't 
opgenomen  fig.  A. 14.26. 

Het  effect  van  het  vrijwel  samenvallen  van  eigenfrequenties  van  "het  dek" 
en  van  "de  vloer"  (gesimplif iceerd  omschreven)  kan  dus  stellig  belangrijk 
zijn  voor  de  doorgeleide  schok  naar  de  vloer. 

Uit  som  8  bleek  verder  dat  er  nog  een  dergelijk  tweetal  trilvormen  aanwezig 
is,  met  iets  verder  uit  elkaar  liggende  eigenfrequenties.  In  de  eerste 
plaats  een  dektrilvorm  met  een  frequentie  van  14.9  Hz  (zie  [4]  fig.  7)  en 
daarnaast  een  trilvorm  van  de  vloer  met  een  frequentie  van  16.5  Hz  (fig. 

8)  . 

De  vraag  rees  of  de  versnellingsniveaux  op  de  vloer  hoger  zouden  zijn 
geweest  als  deze  twee  trilvormen  qua  eigenfrequentie  nog  wat  dichter  bij 
elkaar  hadden  gelegen.  Besloten  werd  om  de  stijfheid  van  het  middelste 
dekveld  zodanig  te  verhogen  dat  de  betreffende  eigenfrequentie  zou  stijgen 
van  14.9  Hz  naar  16.0  Hz.  Hoe  dat  gebeurd  is  en  in  welke  trilvormen  dat 
geresulteerd  heeft  staat  verder  beschreven  in  [4]  onder  som  11. 


Het  eerste  gedeelte  van  het  nu  voorliggende  rapport  betreft  de 
schokresponsie  van  deze  nieuwe  dek  vloer  combinatie  (som  14) . 

Directs  vergelijking  met  de  responsieresultaten  van  som  7  ligt  voor  de  hand 
omdat  het  enige  verschil  in  het  rekenmodel  de  verhoogde  stijfheid  van  het 
middelste  dekveld  betreft. 

Het  tweede  gedeelte  van  dit  rapport  betreft  de  vraag  in  hoeverre  invoering 
van  damping  in  het  model  het  responsiebeeld  gunstig  beinvloedt.  Met  als 
uitgangspunt  de  dek-vloer  combinatie  van  som  7  is  nu  damping  toegevoegd  aan 
het  dek  en  is  de  responsie  opnieuw  bepaald  (som  15) . 


*  *  * 

Van  de  zijde  van  de  K.M. ,  Afdeling  Scheepsbouw,  Bureau  SO&O  werd  het 
onder zoek  begeleid  door  ir.  J.A.A.  Vaders . 
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2.  UITGEVOERDE  ANALYSES. 

(deelhoofdstukken  2.1  t/m  2.13  opzettelijk  blank  gelaten.) 

2.14  SON  14. 

2.14.1  LOCATIE. 

Het  F-dek  t.p.v.  de  commandocentrale .  Daarvan  een  dwarsscheepse  doorsnede 
over  de  voile  breedte  van  het  schip,  ongeveer  halverwege  het  compartiment . 


2.14.2  DEK. 

Model lering  algemeen: 

Het  dek  met  een  breedte  van  16.8  m  is  gemodelleerd  als  een  continue  balk 
met  BECOC  elementen.  Daarin  zijn  de  stijfheids-  en  massagegevens  verwerkt 
van  een  dekstrook  met  een  breedte  van  1.8  m  in  de  langsscheepse  richting. 
Deze  keuze  is  zowel  afgestemd  op  de  onderlinge  afstand  van  de 
dwarsversti jvers  onder  het  dek  (1.8  m)  als  op  de  afmetingen  van  de  vloer 
(rooster  met  een  steek  van  0.6  m) . 

Een  doorsnede  van  deze  continue  balk,  met  daarboven  de  vloer  is  geschetst 
in  fig.  14.1.  Daarin  zijn  ook  mogelijke  plaatsen  van  veren  aangegeven. 
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fiff-  14.1  Dwarsdoorsnede  over  dek  en  vloer. 


Aan  weerszijden  van  hart  schip  is  het  F-dek  ondersteund  door  twee 
langsscheepse  schotten  (onderlinge  afstand  7.2  m).In  het  model  is  de 
dekbalk  op  die  plaatsen  opgelegd  verondersteld. 

De  uiteinden  van  de  balk  worden  stijf  verbonden  geacht  met  de  op  zich 
flexibele  scheepswand,  die  in  het  model  doorloopt  tot  het  G-dek  (3  m 
lager) .Hiermee  wordt  als  het  ware  aan  de  uiteinden  van  de  balk  een 
rotatiesti j fheid  toegevoegd. 
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Sti j  fheidsgegevens : 

Staal,  E  =  210  GN/m^  en  V  =  0.3 
Plaatdikte  dek  7  mm. 

Dwarsverstijver  (webspant) :  225  X  7  +  100  X  10  mm 

De  voor  het  traagheidsmoment  meewerkende  plaatbreedte  bedraagt  volgens 
Lloyds  Register  1.038  m  voor  de  twee  korte  velden  en  1.360  m  voor  het 
centrale  middelste  veld. 

De  traagheidsmomenten  worden  daarmee  resp.  6.5  10"^  m^  en  6.8  10'®  m^. 

De  weerstandsmomenten  zijn  resp.  3.29  10'^  m^  en  3.33  10‘'*  m^. 

Het  af schuifoppervlak  bedraagt  1.575  10"^  m^. 

Het  traagheidsmoment  van  de  twee  3  m  lange  balkelementen,  die  de 
scheepswanden  tot  het  G-dek  vertegenwoordigen  is  geschat  op  6.0  10"^  m^. 

Traagheidsgegevens : 

Aannemende  dat  de  plaatbreedte  1.8  m  is  vinden  we  een  doorsnede  oppervlak 
van  15.175  10'^  m^.  Met  een  dichtheid  van  7800  kg/m^  wordt  de  massa  van  het 
balkprofiel  118.3  kg  per  strekkende  meter.  Deze  massa  is  verder  nog  met  20% 
vermeerderd  i.v.m.  gewicht  kabelbanen,  deel  van  het  gewicht  van  de  veren, 
enz .  De  massa  van  de  langsscheepse  spanten  (om  de  0.6  m)  is  niet  in  de 
berekeningen  betrokken,  omdat  ook  de  bijdrage  aan  de  stijfheid  bij  dit  2- 
dimensionale  model  niet  in  rekening  gebracht  is. 

Een  en  ander  resulteert,  bij  knooppunten  om  de  0.6  m,  in  puntmassa's  van 
85.2  kg.  Er  ontstaat  een  diagonaalmatrix  van  massatermen  die  alleen 
gekoppeld  zijn  aan  translaties  in  verticale  richting. 

Tot  zover  is  de  beschrijving  identiek  aan  die  van  de  in  som  7  doorgerekende 
combinatie.  Zoals  in  de  inleiding  al  uiteengezet  is  nu  t.b.v  som  14  de 
stijfheid  van  het  middelste  dekveld  nog  verder  verhoogd  en  wel  in  de  zelfde 
mate  als  bij  som  11  [4] .  Dus  is  het  traagheidsmoment  van  het  middelste  veld 
verhoogd  van  6.8  10'®  m^  tot  9.0  10'®  m^  en  het  af  schuifoppervlak  van  1.575 
10"®  m^  tot  2. -0846  10"®  m®.  Het  Weerstandsmoment  tegen  bulging  is  evenredig 
mee  verhoogd  van  3.33  10"^  m®  tot  4.41  10"^  m®. 

2.14.3  AANSTOTING . 

De  dekbalk  ondergaat  op  4  plaatsen  identieke  voorgeschreven  bewegingen  in 
verticale  richting:  ter  plaatse  van  de  twee  scheepswanden  en  ter  plaatse 
van  de  schotten. 

De  aanstoting  is  met  de  kickoff  puls  zoals  die  bepaald  werd  met  het 
programme  SWISS  [6]  en  zoals  die,  met  precies  de  zelfde  numerieke  waarden 
eerder  werd  voorgeschreven  bij  som  7  (zie  [3]  hoofdstuk  2.7.3  voor  een 
omschri jving) . 

2.14.4  VLOER. 


Vergeleken  met  som  7  [3]  is  de  vloer  en  de  modellering  daarvan  precies  het 
zelfde  gebleven.  De  beschrijving  ervan  is  hier  nog  een  keer  opgenomen. 

Modellering  algemeen: 

De  vloer  met  een  lengte  van  13.2  m  is  gemodelleerd  als  een  continue  balk 
met  BECOC  elementen.  Daarin  zijn  ook  weer  de  stijfheids-  en  massagegevens 
verwerkt  van  een  strook  met  een  breedte  van  1.8m. 

Deze  gegevens  hebben  betrekking  op  de  combinatie  van  de  vloer  en  de  daarop 
geplaatste  kabinetten. 

De  vloer  ligt  symmetrisch  t.o.v.  hart  schip. 

Sti j  fheidsgegevens : 

Staal,  E  =  210  GN/m®  en  V  =  0.3 

Hit  fig.  14.1  blijkt  dat,  binnen  het  model,  sprake  is  van  vier 
af zonderli jke  doorgaande  profielen. 

Dit  zijn  koudgewalste  UNP  80*50*5  profielen,  conform  de  door  de  KMS  ter 
info  afgegeven  tokening  "voorstel  schokvloeren"  d.d.  31  Mei  1996.  Deze 
tokening  betreft  de  "computerruimte  Sewaco  achter"  Op  verzoek  van  de 
Afdeling  Scheepsbouw  DMKM  is  aangenomen  dat  deze  profielen  ook  voor  de 
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vloer  in  de  commandocentrale  in  aanmerking  komen. 

Van  de  twee  buitenste  profielen  in  fig  14.1  mag  slechts  de  helft  van  de 
buigstijfheid  toegekend  worden  aan  het  balkmodel. 

In  feite  betreft  het  dus  3  parallel  geplaatste  profielen  met  een 
gezamenlijk  traagheidsmoment  van  229.5  10'®  m^  en  een  afschuifoppervlak  van 
1.2  10'^  m^. 

Vergeleken  met  de  aluminivim  palletvloer  in  de  sommen  3,  4  en  6  is  de 
buigstijfheid  El  met  ongeveer  een  factor  7  afgenomen. 

Deze  sti j fheidsgegevens  zijn  niet  van  toepassing  over  de  gehele  lengte  van 
de  vloer,  omdat  er  lokaal  sprake  is  van  een  aanzienlijk  grotere  stijfheid 
veroorzaakt  door  lokaal  aanwezige  kabinetten  die  vast  gemonteerd  zijn  op  de 
vloer.  Uit  fig.  14.2  blijkt  dat  er  sprake  is  van  3  kabinetten,  elk  met  een 
lengte  van  3.6  m  (6  velden)  en  een  breedte  van  1.8  m  (3  velden) .Die 
kabinetten  zijn  symmetrisch  geplaatst  t.o.v  hart  schip  met  tussenruimten 
van  1.2  m  (2  velden) . 


1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  18  10  20  21  22  23 

VLOER  ALS  ASKA  SUBSTR.  2 


KABINET  1 


KABINET  2 


KABINET3(3delen) 


5  7  9  11  13  15  17  19  21  23  25  27  29  31  33  35  37  39  41  43  45  47  49 


1  2  3  4  6  8  10  12  t  P  16  18  20  22  24  26  28  30  32  34  36  t  P40  42  44  46  48  50  51  52 


HUID 


SCHOT 

SHOCK3D-MODEL 


SCHOT  HUID 


?•••••••••  •••••  •-» 

2  3  45  6  7  8  8  10  11  12  13  14  15  16  17  18  19  20  21  22  23  24  25  26  27  28  29  3D 

•  1  DEK  ALS  ASKA  SUBSTR.  1  3i  • 


fig.  14.2  Model  met  knooppuntsnummers  in  SHOCK3D  en  ASKA. 


De  kabinetten  1  en  2  bevatten  elk  doorgaande  profielen  in  langsrichting. 
Over  de  kabinetten  zelf  is  qua  constructiegegevens  nog  niet  voldoende 
bekend  om  die  additionele  stijfheid  met  enige  zekerheid  in  rekening  te 
kunnen  brengen.  Aangenomen  is  dat  aan  de  uiteinden  van  het  kabinet  (het 
eerste  en  het  zesde  veld)  geen  sprake  is  van  een  bijdrage  aan  de  stijfheid 
maar  dat  over  de  vier  middelste  velden  het  traagheidsmoment  en  het 
afschuifoppervlak  een  factor  10  hoger  zijn. 

Het  derde  kabinet  bestaat  in  feite  in  langsrichting  uit  3  af zonderli jke 
delen,  elk  met  een  lengte  van  1.2  m  (2  velden) .  Aangenomen  is  dan  ook  dat 
dit  derde  kabinet  geen  enkele  bijdrage  aan  de  stijfheid  van  de  vloer 
levert . 

Dit  derde  kabinet  brengt  dus  asymmetrie  in  het  model,  niet  alleen  vanwege 
de  stijfheid  van  de  vloer,  maar  ook  vanwege  de  verdeling  van  de  veren, 
zoals  verderop  blijken  zal. 

De  UNP  profielen  vormen  een  volkomen  regelmatig  gelast  frame  met  velden  van 
0.6m  *  0.6m. 


TNO-Rapport 

96-CMC-R1995 


Datum 

16  december  1996 


Pagina 

10 


Traagheidsgegevens : 

Met  een  dichtheid  van  staal  van  7800  kg/m^  is  aan  de  hand  van  de 
beschikbare  afmetingen  geconcludeerd  dat  drie  velden  naast  elkaar  (0.6  m 
bij  1.8  m  )  een  massa  van  24  kg  vormen.  Dat  is  vrijwel  de  zelfde  massa 
zoals  die  eerder  voor  de  palletvloer  [1]  [2]  in  rekening  was  gebracht.  Dat 
was  25  kg.  Die  laatste  waarde  is  hier  gehandhaafd.  Dus  is  aan  elk  van  de 
knooppunten  in  het  model  (onderlinge  afstand  0.6  m)  deze  massa  toegekend, 
met  uitzondering  van  de  knooppunten  aan  de  uiteinden  waar  12.5  kg  is 
gedacht . 

Vender  weegt  elk  van  de  drie  kabinetten  1020  kg.  Die  massa  is  steeds  over  7 
knooppunten  verdeeld,  resp.  85.5  kg  ,  5  maal  170  kg  en  85.5  kg. 

Eventuele  verdere  massa, s  zijn  niet  in  rekening  gebracht. 


2.14.5  VEREN. 

Omdat  ook  aan  de  verende  opstelling  niets  is  veranderd  wordt  hier  voor 
volledige  informatie  verwezen  naar  de  beschrijving  onder  som  7  {  [3]  , 
hoofdstuk  2.7.5). 

Kort  samengevat  is  hier  sprake  van  een  opstelling  op  kabelgewonden  veren 
van  Sebert,  type  K19  100-02  [7]  ,waai~van  de  niet  lineaire  eigenschappen  en 
de  demping  werden  vastgelegd  met  veertype  1  van  het  SHOCK3D  rekenprogramma 
[8]  . 

In  het  onderste  deel  van  fig.  14.3  zijn  in  een  bovenaanzicht  de  plaatsen 
van  de  veren  langs  de  omtrek  van  de  kabinetten  aangegeven.  Per 
dwarsdoorsnede  is  het  aantal  veren  0,  2,  3  of  4.  De  eigenschappen  van  2,  3 
of  4  parallel  gemonteerde  veren  zijn  bij  de  modelvorming  gecombineerd  tot 
enkele  veren  met  de  nummers  l,resp.  2,  resp.  3.  (ook  aangeduid  boven  in 
fig.  14.3) . 


i 

^  Verloop  van  het 

traagheidsmoment 
van  de  vloerbalk 

3  1  1  1  1  1  3  3  ]  1  1  il  1  3  f  ^  ^  1  i  1  2  - 

Veernummers 


1  2  3  4  5  6  7  8  9  10  11  12  13  14  IS  16  17  18  19  20  21  - 

nummers  v.d. 
nietl.  el. 

4222224  4222224  3232323  ^ _ 

aantallen  veren 


^  plaats  v.d.  veren 


fig.  14.3  De  vloerbalk  en  de  daaraan  verbonden  veren. 


Dus  op  21  plaatsen  langs  de  vloerbalk  zijn  in  het  model  niet-lineaire 
elementen  aangebracht . 
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Als  we  de  vergaande  simplif icatie  doorvoeren  dat  onder  "de" 
opstellingsfrequentie  wordt  verstaan  de  freguentie  waarmee  de  totale  massa 
van  de  vloer  en  de  daarop  opgestelde  kabinetten  beweegt  op,  in  dit  geval  54 
parallel  geplaatste  veren,  dan  komt  de  opstellingsfrequentie  uit  op  4.1  Hz. 
De  gemiddelde  statische  veerbelasting  van  een  enkele  veer  is  hier  67  kg, 
hetgeen  binnen  het  opgegeven  [9]  gebied  ligt  van  52  -  104  kg. 


2.14.6  BEREKENINGEN . 

Achtereenvolgens  werden  de  volgende  berekeningen  uitgevoerd: 

(1)  Berekening  van  de  sti j fheidsmatrices  van  het  dek  en  van  de  vloer  als 
aparte  substructures  in  ASKA.  Statische  condensatie  van  het  aantal 
vri jheidsgraden  waardoor  alleen  de  verticale  translaties  overgebleven 
zi  jn. 

(2)  Responsieberekening  van  het  totale  systeem  in  SHOCK3D  [8]  over  500  ms 
met  5000  tijdstappen  van  0.1  ms. 

Uitvoer  van  de  responsie  om  de  10  ms . 

Uitvoer  naar  files  elke  ms  voor  het  plotten  (m.b.v.  TECPLOT)  van 
responsie  resultaten,  voor  het  berekenen  van  de  snedekrachten  en  de 
rotaties  in  de  hierop  volgende  ASKA  berekeningen  en  voor  het 
berekenen  van  schokspectra . 

(3)  ASKA  berekeningen  voor  500  belastingsgevallen  (tijden) . 

Uitvoer  van  de  extreme  waarden  van  de  snedekrachten  in  alle 
knooppunten  met  de  bijbehorende  tijdstippen. 

Uitvoer  van  alle  snedekrachten,  translaties  , rotaties  en 
reaktiekrachten  t.p.v.  de  voorgeschreven  vri jheidsgraden  op  een 
aantal  geselecteerde  tijdstippen  waarop  sprake  is  van  maximale 
responsies .  M.b.v.  IDEAS  zijn  voor  deze  tijdstippen  een  aantal 
verplaatsingsplots  getekend. 

(4)  Berekening  van  initiele  schokspectra  m.b.v.  MATLAB,  op  een  wijze 
identiek  aan  die  bij  Som  7  (zie  figuren  A. 14. 27  -  A. 14. 35). 


2.14.7  REKENRESULTATEN . 

Bij  de  bespreking  van  de  rekenresultaten  zal  veelvuldig  verwezen  worden 
naar  de  knooppuntsnummers  (n.p.  of  node  )  zoals  die  in  ASKA  dan  wel  in 
SH0CK3D  toegekend  werden.  Zie  daartoe  fig.  14.2. 

De  nummers  van  de  niet-lineaire  elementen  (1  t/m  21)  vindt  men  terug  in 
fig.  14.3. 

Alle  figuren  die  betrekking  hebben  op  de  responsie  zijn  in  Appendix  A 
opgenomen . 

2.14.7.1  Schokbeweging  van  het  dek . 

De  absolute  verplaatsingen  van  SHOCK3D  n.p  1  (scheepswand)  en  van  n.p.  6  en 
26  (middens  van  de  velden)  zijn  afgebeeld  in  de  figuren  A. 14.1  en  A. 14. 2  en 
de  versnellingen  van  deze  zelfde  punten  in  de  figuren  A. 14. 3  en  A. 14. 4. 

Deze  4  figuren  betreffen  de  responsie  over  500  ms. 

De  snelheden  van  deze  zelfde  punten  1,  6  en  26  zijn  bovendien  voor  de 
eerste  100  ms  weergegeven  in  fig  A. 14. 5. 

Commentaar : 

Ten  opzichte  van  som  7  is  het  middendek  stijver  maar  de  aanstoting  van  het 
dek  is  niet  veranderd  en  ook  is  "de"  opstellingsfrequentie  van  de  vloer 
niet  veranderd. 

De  schokresponsie  van  het  dek  verschilt  weinig  met  die  van  som  7  alleen  de 
frequentie  van  de  trilling  is  zoals  verwacht  lets  toegenomen  van  circa  33 
Hz  tot  ca  34  Hz  voor  de  zijdekken  en  van  ca  15  Hz  tot  16  Hz  in  het 
middendek.  Dat  blijkt  als  men  bovengenoemde  figuren  vergelijkt  met  de 
figuren  A. 7.1  t/m  A. 7. 5  [3] . 

Die  zelfde  lets  hogere  frequenties  vindt  men  ook  terug  in  som  11  [4] ,  waar 
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16.0  Hz  en  33.8  Hz  gevonden  werd. 

Het  hier  volgende  commentaar  bij  de  dekresponsie  is  dan  ook  vrijwel 
hetzelfde  als  bij  som  7. 

Bij  het  dekveld  tussen  de  scheepshuid  en  het  langsschot  domineert  een 
trilvorm  met  een  eigenfrequentie  van  34  Hz.  De  maximale  doorbuiging  van  dit 
stuk  dek  (n.p.  6)  is  ongeveer  14  mm. 

Het  dek  tussen  de  twee  langsschotten  in  (n.p.  26)  buigt  ongeveer  19  mm 
door.  Dit  is  minder  dan  in  Som  7  (was  23  mm)  door  de  grotere  stijfheid  van 
het  middendek.  Van  de  hier  aanwezige  trilvormen  is,  kijkend  naar  de 
verplaatsingen,  vooral  die  met  een  frequentie  van  16  Hz  belangrijk. 

De  maximum  waarden  van  de  versnellingen  op  het  dek  (fig  A. 14. 3  en 
fig. A. 14. 4)  liggen  vrijwel  steeds  (aanzienlijk)  beneden  de  1000  m/s^. 
Blijkbaar  worden  ook  hogere  trilvormen  in  het  dek  aangestoten,  hetgeen  ook 
verwacht  mag  worden  gezien  de  kleine  stijgtijd  tot  de  maximum  snelheid  van 
de  voorgeschreven  puls . 

Uit  de  ti jdsignalen,  met  name  de  versnellingen  blijkt  dat  het  systeem  licht 
gedempt  is.  De  amplituden  nemen  langzaam  af .  Damping  is  alleen  in  het  model 
ingevoerd  in  de  kabelgewonden  veren.  Door  allerlei  andere  oorzaken  zal  de 
damping  in  werkeli jkheid  groter  zijn  en  zullen  met  name  de  trillingen  met 
hogere  frequenties  sneller  uitgedempt  zijn. 

2.14.7.2  Schokbeweging  van  de  vloer. 

De  versnellingen  van  SHOCK3D  n.p.  5,  11,  17,  27,  39,  43  en  47  zijn 
afgebeeld  in  de  figuren  A. 14. 6  t/m  A. 14. 12.  Bovendien  zijn  van  de 
knooppunten  7,  11  en  15  (onder  kabinet  1)  de  versnellingen  weergegeven  in 
fig.  A. 14. 13  en  van  de  knooppunten  21,  27  en  33  (onder  kabinet  2)  in  fig. 

A.  14. 14  en  van  de  knooppunten  37,  43  en  49  (onder  het  gedeelde  kabinet  3  ) 
in  fig.  A. 14. 15.  Helaas  was  in  [3]  de  figuur  met  versnellingen  van  de 
knooppunten  37,  43  en  49  weggevallen.  Deze  figuur  wordt  alsnog  gegeven  in 
figuur  A. 14. 26. 

Verder  zijn  op  een  zestal  kritieke  tijdstippen  de  vloerverplaatsingen 
weergegeven.  Zie  de  figuren  A. 14. 16  t/m  A. 14. 21. 

Commentaar : 

Vergeleken  met  som  7  is  het  middendek  stijver  geworden. 

De  vloer  laat  nu  een  wat  andere  schokresponsie  zien. 

Alvorens  daar  op  in  te  gaan  kan  opgemerkt  worden  dat  in  alle 
versnellingssignalen  opnieuw  de  laagfrekwente  opstellingsfrequentie  van  4.1 
Hz  herkenbaar  is.  Die  laagfrekwente  component  is  verantwoordeli jk  voor  een 
versnellingsniveau  over  de  gehele  lengte  op  de  vloer  van  ongeveer  30  m/s^. 
Daarnaast  zijn  versnellingscomponenten  met  hogere  frequenties  aanwezig. 

Met  namet.p.v.  de  buitenste  knooppunten  5  en  17  van  kabinet  1  dragen  deze 
versnellingen  weer  aanzienlijk  bij  aan  het  totale  versnellingsniveau  tot 
ongeveer  90  resp.  70  m/s^.  Dit  is  dus  iets  lager  dan  bij  som  7. 

Een  frequentie  van  ongeveer  33  Hz,  dicht  bij  die  van  34  Hz  geconstateerd  in 
de  buitenste  dekvelden,  blijkt  hier  in  sterke  mate  opgebouwd  te  worden. 

Als  we  nu  kijken  naar  de  versnellingssignalen  van  de  knooppunten  5,  7,  11, 
15  en  17  onder  dit  kabinet  1  dan  is  ook  nu  het  signaal  van  het  middelste 
knooppunt  11  in  tegenfase  met  de  signalen  van  de  vier  andere  knooppunten. 
Blijkbaar  treedt  er  ook  nu  lokaal  een  buigtrilling  in  de  vloer  onder 
kabinet  1  op.  Dat  blijkt  natuurlijk  ook  uit  de  berekende  trilvormen  [4]  van 
de  vloer.  Er  is  nog  steeds  sprake  van  een  eigentrillingsvorm  van  de  vloer 
([4],  fig  33)  die  met  een  eigenfrequentie  van  31.7  Hz  dicht  ligt  bij  de 
eigenfrequentie  van  33.8  Hz  van  een  dektrilling  als  excitatiebron.  Hierdoor 
ontstaat  nu  een  zwevingsfrequentie  van  2.1  Hz  bovenop  de  gemiddelde 
frequentie  van  33  Hz.  Daardoor  is  de  trilling  na  0.25  s  maximaal  en  lijkt 
ze  na  0.5  s  uitgedempt  te  zijn.  In  som  7  waren  deze  eigenfrequenties  32.7 
resp.  31.7  Hz.  Hier  ontstond  dus  een  zwevingsfrequentie  van  1  Hz,  die  dus 
na  0.5  s  juist  op  zijn  maximum  is. 

Dat  de  knooppunten  5  en  17  veel  grotere  amplituden  te  zien  geven  dan  de 
tussenliggende  knooppunten  is  niet  zo  verwonderli jk  gezien  de  geringe 
stijfheid  van  de  vloer  tussen  de  knooppunten  5  en  7  en  tussen  de 
knooppunten  15  en  17.  Die  uiteinden  van  het  kabinet  worden  gemakkelijk 
aangestoten  omdat  het  aantal  veren  daar  juist  twee  keer  zo  hoog  is  en  de 
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massa  twee  keer  zo  laag  is,  vergeleken  met  de  tussenliggende  punten.  Een  en 
ander  blijkt  uit  de  vorm  van  de  vloertrilling  ([4],  figuur  33). 

Kabinet  2  is  identiek  aan  kabinet  1,  evenals  de  onderliggende 
vloerconstructie .  Toch  is  hier  de  responsie  van  de  uitgekozen  knooppunten 
21,  27  en  33  met  piekversnell ingen  tot  ongeveer  40  m/s^  veel  gunstiger  dan 
bij  kabinet  1.  Die  knooppunten  liggen  echter  toevallig  alle  drie  vrijwel 
precies  in  de  knopen  van  die  31.7  Hz  trilvorm. 

De  tussenliggende  knooppunten  23  ,25,  29  en  31  zullen  wel  een  bijdrage  te 
zien  geven  van  die  trilvorm. 

Voor  wat  betreft  de  onder  kabinet  2  geselecteerde  en  gepresenteerde 
signalen  is  het  verschil  met  som  7  erg  klein. 


Het  derde  kabinet  werd  opgebouwd  gedacht  uit  3  af zonderli jke  losstaande 
delen,  dus  zonder  bij  te  dragen  aan  de  stijfheid  van  de  vloer. 

In  dit  rechter  deel  van  de  vloer  zal  de  eerder  genoemde  31.7  Hz  trilvorm 
van  de  vloer  dus  een  geheel  andere  vorm  hebben  dan  onder  kabinet  1 . 

Dat  blijkt  ook  wel  uit  [4],  fig. 33.  De  amplituden  zijn  wat  minder  groot  dan 
onder  het  kabinet  1  en  dat  is  ook  het  beeld  wat  uit  de  berekende  responsie 
naar  voren  komt.  Ook  hier  blijkt  (fig  A. 14. 15  )  een  sterke  33  Hz  component 
aanwezig.  De  maximum  versnellingspieken  lopen  hier  op  tot  bijna  70  m/s^. 
Ook  hier  valt  de  zwevingsfrequentie  van  ongeveer  2  Hz  op  waardoor  de 
responsie  van  de  vloer  na  0.5  s  klein  is. 

De  grootste  versnelling  die  op  de  vloer  optreedt  bedraagt  dus  ongeveer 
90  m/s^,  waarvan  ongeveer  60  m/s^  gerelateerd  is  aan  een  trilling  met  een 
frequentie  van  ongeveer  33  Hz.  Die  laatste  trilling  gaat  dan  gepaard  met 
een  verplaatsingsamplitude  van  1.3  mm. 

Tenslotte  wordt  geconcludeerd  dat  de  responsie  van  de  vloer  maar  in  zeer 
geringe  mate  blijk  geeft  aangestoten  te  worden  in  de  16.5  Hz  trilvorm  (zie 
[4],  fig.  30)  of  in  de  daar  vlakbij  gelegen  trilvorm  met  een  frequentie  van 
16.0  Hz  (zie  [4],  fig.  29).  Dat  is  ook  niet  zo  verwonderli jk  omdat  bij 
beide  trilyormen  het  middelste  dekveld  in  tegenfase  met  de  twee  zijvelden 
beweegt.  Dit  in  tegenstelling  tot  het  eerder  genoemde  en  dominante  paar 
trilvormen  (31.7  en  33.8  Hz)  waarbij  alle  drie  de  dekvelden  in  fase 
bewegen.  Een  zweving  in  de  responsie  tussen  de  16.0  Hz  en  de  16.5  Hz 
trilvormen  kon  dan  ook  niet  worden  gesignaleerd. 


2.14.7.3  Gedrag  van  de  veren . 

Voor  de  niet-lineaire  elementen  2,  4,  6,  en  11  (zie  fig.  14.3  voor  de 
nummering)  zijn  de  gedurende  500  ms  doorlopen  veerkarakteristieken 
weergegeven  in  de  figuren  A. 14. 22  t/m  A. 14. 25. 

Commentaar : 

Elk  van  de  4  elementen  vertegenwoordigt  2  parallel  geplaatste  K19  100-02 
veren . 

De  elementen  2,  4  en  6  zijn  alle  onder  kabinet  1  geplaatst. 

De  veren  blijken  maximaal  ongeveer  40  mm  ingedrukt  en  ook  uitgerekt  te 
worden.  Ter  plaatse  van  de  langsschotten  zijn  deze  waarden  het  kleinst. 
Voor  deze  niet  bij zonder  grote  waarden  mag  worden  aangenomen  dat  de 
veerkarakteristiek  qua  niet-lineariteit  goed  in  het  model  verwerkt  is.  De 
verschillen  met  som  7  blijken  klein. 

2.14.7.4  Buigspanningen  in  het  dek  en  in  de  vloer. 

Voor  500  quasi-statische  belastingsgevallen  is  m.b.v.  ASKA  per  knooppunt 
zowel  voor  het  dek  als  voor  de  vloer  nagegaan  wat  de  extreme  waarden  zijn 
voor  de  buigende  momenten  in  de  balkelementen .  Tabel  2.3.1  geeft  daarvan 
een  overzicht  met  de  grootste  positieve  en  de  grootste  negatieve  waarden 
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evenals  de  tijdstippen  waarop  die  extremen  optreden. 

Omdat  bij  het  dek  duidelijk  sprake  is  van  symmetrie  (de  extreme  waarden 
meestal  tot  op  3  cijfers  gelijk  zijnde)  vermeldt  de  tabel  alleen  de  waarden 
voor  de  linker  helft  van  de  constructie. 

V.w.b.  de  vloer  zijn  in  de  tabel  niet  die  knooppunten  opgenomen  die 
verbonden  zijn  aan  de  verstijfde  balkelementen  { samenhangend  met  de 
aanwezigheid  van  kabinetten) .  Daar  is  de  situatie  met  dit  voorlopige  model 
minder  goed  interpreteerbaar . 

Als  we  de  buigende  momenten  delen  door  het  weerstandsmoment  vinden  we  dat 
in  het  dek  de  grootste  buigspanning  optreedt  na  112  ms  in  knooppunt  10  (bij 
het  langsschot)  met  495  MPa,  aanzienlijk  hoger  dan  in  som  7  (402  MPa) .  De 
op  een  na  grootste  spanning  treedt  op  in  knooppunt  6  (hart  linker  dek)  na 
140  ms  en  bedraagt  326  Mpa.  De  spanningen  in  het  middendek  zijn  door  het 
hogere  buigweerstandsmoment  aanzienlijk  gedaald  (maximaal  263  MPa) . 

De  vloeigrens  van  het  toegepaste  staal  ligt  met  ongeveer  350  MPa  behoorlijk 
lager  dan  de  hoogst  gevonden  waarde  en  een  elastische  berekening,  zoals 
hier  uitgevoerd,  wordt  dan  ook  discutabel. 

In  de  UNP  profielen  van  de  vloer  treedt  een  grootste  buigspanning  op  van  87 
MPa  in  knooppunt  14  na  289  ms.  Ongeveer  10%  minder  grote  spanningen  treden 
op  in  de  knooppunten  2  en  10.  Vergeleken  met  de  vloeigrens  is  het 
spanningsniveau  in  de  vloer  laag  en  vergeleken  met  som  7  ca  25%  lager.  Ook 
bij  deze  berekening  treden  de  grootste  spanningen  in  de  vloerbalk  weer  op 
die  plaatsen  op  waar  sprake  is  van  een  overgang  in  buigsti j fheid. 
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Extreme 

waarden 

Positief 

Negatief 

Mb  X  10^ 

t 

Mb  X  10® 

t 

[Nm] 

[ms] 

[Nm] 

[ms] 

0.74 

38 

0.83 

23 

0.29 

42 

0.29 

74 

0.79 

23 

0.71 

38 

0.94 

21 

0.82 

37 

1.04 

140 

0.97 

41 

0.91 

85 

0.72 

41 

0.60 

18 

0.47 

44 

0.66 

12 

0.70 

115 

1.43 

69 

1.63 

112 

0.97 

9 

0.91 

52 

0.64 

9 

0.63 

203 

0.53 

47 

0.53 

14 

0.70 

114 

0.55 

9 

1.00 

109 

0.86 

9 

1.16 

108 

1.01 

8 

0.046 

135 

0.045 

242 

0.039 

139 

0.020 

215 

0.025 

273 

0.021 

350 

0.014 

241 

0.017 

320 

0.015 

228 

0.013 

334 

0.033 

412 

0.038 

304 

0.039 

210 

0.050 

289 

0.022 

210 

0.017 

136 

0.020 

227 

0.016 

148 

0.027 

227 

0.021 

149 

0.011 

404 

0.012 

71 

0.034 

149 

0.031 

227 

0.030 

148 

0.028 

227 

0.037 

135 

0.032 

213 

0.028 

133 

0.026 

212 

Knooppunt- 

nuitimers 

in 

ASKA 


Weerstands 

moment 

Wb  X  10'" 


vloer 


0.574 


Tabel  2.14.1 


De  extreme  buigende  momenten  (som  14) 
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2.15  SOM  15. 


2.15.1  Geqevens . 

Som  15  is  wat  betreft  lokatie,  samenstelling  dek  en  vloer,  aanstoting, 
veren  tussen  dek  en  vloer  gelijk  aan  Som  7  en  het  verschil  met  de  vorige 
Som  14  is  slechts  het  feit  dat  hier  weer  de  oorspronkeli jke  stijfheid  van 
het  middendek  is  gebruikt.  De  bedoeling  van  de  huidige  som  is  het 
onderzoeken  van  het  effect  van  damping  op  de  opslingering  in  versnellingen 
die  in  Som  7  was  geconstateerd. 

Damping  in  het  dek  is  geintroduceerd  door  tussen  alle  knooppunten  van  het 
dek  horizontaal  liggende  veren  te  plaatsen,  zonder  stijfheid,  maar  met 
alleen  een  waarde  voor  de  demping  in  de  verticale  richting.  Deze  methode  is 
eerst  getest  op  eenvoudiger  opstellingen  (een  dekveld  en  alleen  het  dek) . 
Zie  Appendix  B  voor  een  meer  uitgebreide  beschrijving.  In  deze  som  zijn 
dezelfde  veren  gebruikt  als  in  deze  voorbereidende  sommen.  De 
dempingsfactor  van  de  veren  van  veertype  1  in  het  SH0CK3D  rekenprogramma 
[8]  is  K=6500  Ns/m.  Deze  dempingsfactor  geeft  voor  de  belangri jkste 
trilvormen  een  demping  van  3  a  4%  van  de  kritische  demping. 

2.15.2  BEREKENINGEN. 

Achtereenvolgens  warden  de  volgende  berekeningen  uitgevoerd: 

(1)  Berekening  van  de  sti j fheidsmatrices  van  het  dek  en  van  de  vloer  als 
aparte  substructures  in  ASKA.  Statische  condensatie  van  het  aantal 
vri jheidsgraden  waardoor  alleen  de  verticale  translaties  overgebleven 
zijn. 

(2)  Responsieberekening  van  het  totale  systeem  in  SH0CK3D  [11]  over  500 
ms  met  5000  tijdstappen  van  0.1  ms. 

Uitvoer  van  de  responsie  om  de  10  ms. 

Uitvoer  naar  files  elke  ms  voor  het  plotten  (m.b.v.  TECPLOT)  van 
responsie  resultaten  en  voor  schokspectra . 

(3)  Berekening  van  initiele  schokspectra  m.b.v.  MATLAB,  op  een  wijze 
identiek  aan  die  bij  Som  7  [4]  en  Som  14  (zie  figuren  A.  15. 20  - 
A. 15. 27)  . 

2.15.3  REKENRESULTATEN . 

Bij  de  bespreking  van  de  rekenresultaten  wordt  weer  veelvuldig  verwezen 
naar  de  knooppuntsniimmers  zoals  die  in  ASKA  dan  wel  in  SHOCK3D  toegekend 
werden.  Zie  daartoe  fig.  14.2. 

De  nummers  van  de  niet-lineaire  elementen  (1  t/m  21)  vindt  men  terug  in 
fig.  14.3. 

Alle  figuren  die  betrekking  hebben  op  de  responsie  zijn  in  appendix  A 
opgenomen . 

2.15.3.1  Schokbeweging  van  het  dek . 

De  absolute  verplaatsingen  van  SHOCK3D  n.p  1  (scheepswand)  en  van  n.p.  6  en 
26  (middens  van  de  velden)  zijn  afgebeeld  in  de  figuren  A. 15.1  en  A. 15. 2  en 
de  versnellingen  van  deze  zelfde  punten  in  de  figuren  A. 15. 3  en  A.15.4. 

Deze  4  figuren  betreffen  de  responsie  over  500  ms. 

De  snelheden  van  deze  zelfde  punten  1,  6  en  26  zijn  bovendien  voor  de 
eerste  100  ms  weergegeven  in  fig  A. 15. 5. 

Commentaar : 

Ten  opzichte  van  som  7  is  de  invloed  van  de  ingevoerde  demping  zeer 
duidelijk  in  de  zijdekken  (knooppunt  6) .  Terwijl  in  som  7  een  beweging  werd 
gevonden  met  een  voornaamste  frequentie  van  ca  33  Hz  en  vri j  wel  constante 
amplitude  zien  we  nu  dat  de  verplaatsing  na  0.5  s  al  vrijwel  constant  is. 
Het  middendek  blijft  langer  trillen  omdat  lagere  freguenties  minder  snel 
uitdempen.  Wel  is  te  zien  dat  de  hogere  freguenties  sneller  verdwijnen  in 
vergelijking  met  Som  7.  Dit  laatste  is  zeker  goed  te  zien  in  de 
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resulterende  versnell ingen.  De  optredende  maximale  versnelling  is  minder 
dan  in  Som  7  (ca  900  m/s^  in  plaats  van  1200  m/s^)  en  treedt  nu  nog  maar  een 
keer  op.  Met  snelheids  signaal  lijkt  veel  op  dat  van  Som  7  maar  is  door  het 
uitdempen  van  de  hoge  frequenties  veel  rustiger. 

Een  en  ander  blijkt  als  men  bovengenoemde  5  figuren  vergelijkt  met  de 
figuren  A. 7.1  t/m  A. 7. 5  [3] . 

De  maximale  doorbuiging  van  de  zijdekken  van  13  mm  is  iets  kleiner  dan  in 
Som  7  {14  mm) .  Het  middendek  trilt  weer  in  dezelfde  frequentie  van  ca  15  Hz 
met  een  maximale  doorbuiging  van  15.6  mm,  aanzienlijk  minder  dan  in  Som  7 
(23  mm) . 

2.15.3.2  Schokbeweging  van  de  vloer. 

De  versnellingen  van  SHOCK3D  n.p.  5,  11,  17,  27,  39,  43  en  47  zijn 
afgebeeld  in  de  figuren  A. 15. 6  t/m  A. 15. 12.  Bovendien  zijn  van  de 
knooppunten  7,  11  en  15  (onder  kabinet  1)  de  versnellingen  weergegeven  in 
fig.  A. 15. 13  en  van  de  knooppunten  21,  27  en  33  (onder  kabinet  2)  in  fig. 

A. 15. 14  en  van  de  knooppunten  37,  43  en  49  in  fig.  A. 15. 15. 

Commentaar : 

De  invloed  van  de  ingevoerde  dekdemping  is  ook  in  de  beweging  van  de  vloer 
merkbaar.  Zo  zien  we  bij  knooppunt  5  waar  in  Som  7  de  versnelling 
opslingerde  tot  ca  100  m/s^,  dat  nu  nog  slechts  ca  65  m/s^  wordt  bereikt. 
Aangezien  geen  damping  is  aangebracht  in  de  vloer  (het  is  heel  goed 
mogelijk  dat  deze  minder  gedempt  wordt  dan  het  dek)  zien  we  geen  bijzonder 
sterke  afname  in  de  trilling.  Maar  dat  hangt  ook  samen  met  het  optreden  van 
een  zweving  die  bij  de  huidige  SOM  15  weer  een  frequentie  van  1  Hz  heeft, 
zodat  we  na  0.5  sec  juist  op  het  maximum  zitten.  (was  2  Hz  bij  SOM  14) . 

Ook  op  de  andere  knooppunten  zien  we  een  reductie  van  de  maximale 
versnellingen. 

Wederom  zien  we  in  alle  versnellingssignalen  een  laagfrequente 
opstellingsfrequentie  van  4.1  Hz.  Onder  kabinet  1  overheerst  weer  de 
beweging  met  een  frequentie  van  33  Hz,  die  verantwoordeli jk  is  voor  een 
aanzienlijke  verhoging  van  het  versnellingsniveau,  vooral  bij  de  buitenste 
knooppunten  5  (65  m/s^)  en  17  (55  m/s^)  van  kabinet  1. 

Onder  kabinet  2  is  de  situatie  weer  gunstiger  omdat  hier  vanuit  het  dek  de 
lagere  frequentie  van  15  Hz  overheerst.  De  versnellingen  zijn  hier  vrijwel 
gelijk  aan  die  van  Som  7,  alleen  ontbreken  de  hoogfrequente  piekjes. 

Ook  onder  kabinet  3  verschillen  de  versnellingen  niet  veel  met  die  van 
Som  7.  De  extrema  zijn  vrijwel  gelijk  en  de  amplitude  van  de  33  Hz  trilling 
is  iets  kleiner. 

2.15.3.3  Gedrag  van  de  veren . 

Voor  de  niet-lineaire  elementen  2,  4,  6,  en  11  (zie  fig.  14.3  voor  de 
nummering)  zijn  de  gedurende  500  ms  doorlopen  veerkarakteristieken 
weergegeven  in  de  figuren  A. 15. 16  t/m  A. 15. 19. 

Commentaar : 

De  verschillen  met  Som  7  zijn  klein. 
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3 .  CONCLUSIES . 

(1)  Vergeleken  met  een  eerder  doorgerekende  dek-vloer  combinatie  met 
staler.  UNP  profielen  (som  7)  t.p.v.  de  commandocentrale  van  het  LCF 
[3]  is  nu  met  Som  14  het  effect  van  ongeveer  samenvallende 
eigenfrequenties  bekeken  en  met  Som  15  het  effect  van  een  demping  van 
3  a  4%  in  het  dek. 

Vender  geldt: 

De  dekconstructie  is  onveranderd  gelaten,  uitgezonderd  het 
middenveld  dat  nu  bij  Som  14  stijver  is. 

De  aanstoting  van  het  dek  met  de  kickoff  puls  t.p.v.  de  huid  en 
de  langsschotten  is  conform  Som  4  en  Som  7 . 

De  vloer  bestaat  evenals  in  Som  7  uit  een  gelast  staler  frame 
van  UNP  80*50*5  profielen. 

Evenals  in  Som  7  zijn  staalkabelveren  van  Sebert  K19  100-02 
gebruikt  tussen  dek  en  vloer. 

(2)  De  schokresponsie  van  het  dek  is  in  Som  14  door  de  hogere  stijfheid 
van  het  middendek  iets  anders.  Dit  uit  zich  met  name  in  iets  andere 
eigenfrequenties  van  de  twee  belangri jkste  trilvormen.  De  trilvorm 
waarbij  zijdekken  en  middendek  in  tegenfase  bewegen  heeft  een 
frequentie  van  16  i.p.v.  15  Hz  en  de  trilvorm  waarbij  alle  dekken  in 
fase  bewegen  heeft  een  frequentie  van  33.8  Hz  (was  32,7  Hz). 

(3)  De  schokresponsie  van  de  vloer  is  in  Som  14  opvallend  verschillend 
van  Som  7.  Doordat  de  eigenfreqentie  van  33.8  Hz  verder  bij  een 
eigenfrequentie  van  31.7  Hz  van  de  vloer  vandaan  is  gekomen,  ontstaat 
in  het  versnellingsignaal  een  zweving  met  een  frequentie  van  2  Hz.  In 
Som  7  was  deze  verschilfrequentie  1  Hz .  In  Som  14  zien  we  in  de 
versnellingen  van  de  vloer  een  voile  periode  terug  van  deze  zweving 
en  in  Som  7  slechts  een  halve  periode  gedurende  de  500  ms  waarover  de 
responsie  berekend  is .  Hierdoor  leek  het  dat  bij  Som  7  de 
versnellingen  zouden  blijven  opslingeren. 

(4)  Door  de  grotere  stijfheid  van  het  middendek  zijn  de  spanningen  in  dit 
deel  van  het  dek  lager  dan  in  Som  7 .  Daar  staat  echter  tegenover  dat 
daardoor  de  spanningen  in  de  zijdekken  juist  extra  hoog  zijn  geworden 
(ca  25  %  hoger) .  De  spanningen  in  de  vloer  blijken  juist  ca  25%  lager 
te  zijn. 

(5)  In  Som  15  is  de  schokresponsie  van  het  dek  ten  gevolge  van  de 
ingevoerde  demping  aanzienlijk  verschillend  van  die  bij  Som  7.  Zoals 
te  verwachten  is  worden  vooral  de  dektrilvormen  met  hogere 
frequenties  snel  uitgedempt.  Ook  de  belangri jke  trilvorm  met  de 
frequentie  van  32.7  Hz  is  na  0.5  s  al  vrijwel  verdwenen. 

De  schokresponsie  van  de  vloer  verandert  ten  opzichte  van  Som  7 
minder  sterk.  De  vloer  zelf  is  in  de  berekeningen  niet  gedempt.  Toch 
is  de  opslingering  t.g.v.  de  zweving  aanzienlijk  verkleind  en  zijn  de 
amplituden  van  de  vloertrillingen  onder  de  kabinetten  1  en  3  wel  wat 
afgenomen.  Onder  kabinet  1  zelfs  van  100  tot  ongeveer  65  m/s^. 


(6) 
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fig. 


fig. 


Node  nr  1 


.14.3  Versnelling  som  14  dek  n.p.  6  en  n.p.  1  (voorgeschreven) . 

-  Node  nr  1 


.14.4  Versnelling  som  14  dek  n.p.  26  en  n.p.  1  (voorgeschreven) 
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-  Node  nr  11 

t  [m/s*] 


fig.  A.  14. 7  Versnelling  som  14  vloer  n.p.  11. 

-  Node  nr  17 

T  [m/s*] 


fig.  A.  14. 8  Versnelling  som  14  vloer  n.p.  17. 
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-  Node  nr  27 

t  [m/s*] 


fig.  A. 14. 9  Versnelling  som  14  vloer  n.p.  27. 

-  Node  nr  39 

t  [m/s*] 


fig.  A. 14. 10 


Versnelling  som  14  vloer  n.p.  39. 
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-  Node  nr  43 

T  [m/s"] 


problem  14 


fig.  A. 14. 11  Versnelling  som  14  vloer  n.p.  43. 

-  Node  nr  47 

t  [m/s"] 


problem  14 


fig.  A. 14. 12  Versnelling  som  14  vloer  n.p.  47. 
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problem  14 


fig.  A. 14. 13  Versnelling  som  14  vloer  n.p.  7,  11  en  15. 
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fig.  A. 14. 14 


Versnelling  som  14  vloer  n.p.  21,  27  en  33. 
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fig.  A. 14.20  Verplaatsing  vloer  som  14  na  289  ms. 


fig.  A14.32  Verplaatsing  vloer  som  14  na  304  ms. 


De  verplaatsingen  zijn  in  deze  figuren  een  faktor  50  opgeschaald. 
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problem14 

fig.  A.  14.24  Veerkarakteristiek  som  14  element  6. 
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fig.  A. 14.25 
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Veerkarakteristiek  som  14  element  11. 
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Som  14;  Dek;  N.p.  1 


fig.  A. 14. 27 


Initieel  schokspectriam,  som  14,  dek,  n.p.  1 
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Som  14;  Dek;  N.p.  6 


fig.  A. 14. 28 


Initieel  schokspectriim,  som  14,  dek,  n.p.  6 
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Sm  14;  Dek;  N.p.  26 


fig.  A. 14. 29 


Initieel  schokspectrum,  som  14,  dek,  n.p. 26 


fig.  A. 14. 30 


Initieel  schokspectrum,  som  14,  vloer,  n.p.  5 
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fig.  A. 14. 31 


Initieel  schokspectrum,  som  14,  vloer,  n.p.  11 
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fig.  A. 14. 32 


Initieel  schokspectrum,  som  14,  vloer,  n.p.  17 
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Som  14;  Vloer;  N.p.  27 


fig.  A. 14. 33 


Initieel  schokspectrum,  som  14,  vloer,  n.p.  27 
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Som  14;  Vloer;  N.p.  37 


fig.  A. 14. 34 


Initieel  schokspectrum,  som  14,  vloer,  n.p.  37 
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S(M  14;  Vloer;  N.p.  43 


fig.  A. 14. 35 


Initieel  schokspectrum,  som  14,  vloer,  n.p.  43 
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t[m/s] 


fiff-  A. 15. 5  Snelheid  som  15  dek  n.p.  6,  26  en  1  (voorgeschreven) 


-  Node  nr  5 

t  [m/s*] 


problem  15 


fig.  A. 15. 6  Versnelling  som  15  vloer  n.p.  5. 
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-  Node  nr  27 

T  [m/s*] 


fig.  A. 15. 9  Versnelling  som  15  vloer  n.p.  27. 


-  Node  nr  39 

t  [m/s*] 


fig.  A.  15. 10  Versnelling  som  15  vloer  n.p.  39. 
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-  Node  nr  43 

t  [m/s*] 


fig.  A. 15. 11  Versnelling  som  15  vloer  n.p.  43. 


-  Node  nr  47 

t  [m/s*] 


fig.  A. 15. 12  Versnelling  som  15  vloer  n.p.  47. 
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. Node  nr  37 

t  [m/s*]  -  Node  nr  43 


fig.  A. 15. 15  Versnelling  som  15  vloer  n.p.  37,  43  en  49 
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Som  15;  Dek;  N.p.  6 


fig.  A. 15. 20 


Initieel  schokspectrum,  som  15,  dek,  n.p.  6 
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Som  15;  Dek;  N.p.  26 


fig.  A. 15. 21 


Initieel  schokspectrum,  som  15,  dek,  n.p.  26 


Sam  15;  Vloer;  N.p.  5 


fig.  A. 15. 22 


Initieel  schokspectrum,  som  15,  vloer,  n.p.  5 
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Som  15;  Vloer;  N.p.  11 


fig.  A. 15. 23 


Initieel  schokspectriam,  som  15,  vloer,  n.p.  11 
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Som  15;  Vloer;  N.p.  17 


fig.  A. 15. 24 


Initieel  schokspectrum,  som  15,  vloer,  n.p.  17 
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S(M  15;  Vloer;  N.p.  27 


fig.  A. 15. 25 


Initieel  schokspectrum,  som  15,  vloer,  n.p.  27 
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Som  15;  Vloer;  N.p.  37 


fig.  A. 15. 26 


Initieel  schokspectrum,  som  15,  vloer,  n.p.  37 
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Soi  15;  Vloer;  N.p.  43 


fig.  A.  15. 27 


Initieel  schokspectrum,  som  15,  vloer,  n.p.  43 
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Appendix  B  Demping  van  het  dek 

De  modeliering  van  demping  is  getest  met  een  eenvoudig  model  met  slechts 
het  linker  dekveld  met  aan  beide  zijden  de  "kickoff  puls"  zoals  die  in 
SOTi  7  is  gebruikt.  Dit  model  wordt  verder  niet  gerapporteerd.  Vervolgens  is 
als  testprobleem,  som  15a  met  alleen  het  volledige  dek  doorgerekend .  De 
gegevens  zijn  verder  identiek  aan  die  van  som  7.  Slechts  de  vloer  en  de 
veren  zijn  weggelaten.  De  demping  is  gemodelleerd  door  tussen  alle 
opvolgende  knooppunten  van  het  dek  veren  te  plaatsen  met  als 
afschuifs,tij  fheid  in  de  vertikale  richting: 

Force  =  K  x  (verschil  in  vertikale  snelheid  van  de  aanliggende  knooppunten) 
met  K=6500  Ns/m. 

In  de  volgende  figuren  worden  versnelling,  snelheid  en  verplaatsing  van  de 
middens  van  linker  en  middendekveld  getoond,  samen  met  de  voorgeschreven 

•  Die  van  het  rechter  veld  zijn  vanwege  de  symmetrie  gelijk  aan  die 
van  het  linker  veld. 

- -  Node  nr  1 


£ig«  B.l  Versnelling  som  15a  dek  n.p.  6,  26  en  n.p.  1  (voorgeschreven) . 

De  versnelling  van  het  linker  veld  (knooppunt  6)  toont  al  snel  een  mooie 
gedempte  sinusvorm  met  een  freguentie  van  ca  33  Hz  en  een  demping  van  3% 
zoals  is  te  zien  in  figuur  B.l  .  Deze  freguentie  behoort  bij  een 
dektrilvorm  met  6  knooppunten  waarbij  al  de  dekvelden  in  dezelfde  richting 
bewegen.  De  versnelling  van  het  middenveld  toont  een  combinatie  van 
trillingen,  waarvan  die  met  een  freguentie  van  ca  33  Hz  en  van  ca  14  Hz  de 
belangri jkste  zijn. 

De  snelheden  in  alle  plaatvelden  zijn  complex  doordat  meerdere  trilvormen 
aanwezig  zijn,  zie  figuur  B.2.  Ook  hier  zijn  de  freguenties  van  33  en  14  Hz 
weer  het  belangri jkst . 
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De  verplaatsing  van  het  linker  veld  toont  duidelijk  een  combinatie  van 
trilvormen,  zie  figuur  B.3.  De  verplaatsing  van  het  middenveld  wordt 
overheerst  door  een  trilvorm  met  een  f reguentie  van  ca  14  Hz .  Deze 
frequentie  behoort  bij  een  dektrilvorm  met  4  knooppunten  (de 
aanstootpunten) ,  waarbij  het  middendek  de  tegenovergestelde  verplaatsing 
heeft  van  de  zijvelden.  Indien  de  invloed  van  de  trilvorm  met  33  Hz  wordt 
verwaarloosd,  volgt  een  demping  voor  de  hoofdtrilvorm  van  4%.  Er  kan 
opgemerkt  worden  dat  in  het  begin  van  de  beweging  alle  dekken  in  fase 
bewegen  (de  trilvorm  met  33  Hz) ,  terwijl  zij  later  in  tegenfase  gaan 
bewegen  (de  trilvorm  met  14  Hz) . 

Hier  volgt  een  korte  samenvatting  van  de  theorie . 

Voor  een  een-massa-veer  systeem  is  de  bewegingsvergelijking: 

M  X  +  C  X  +  K  X  =  F 


waarin 

x:  verplaatsing  van  de  massa; 

M:  massa; 

C:  demping; 

K;  stijfheid  van  de  veer; 

F:  uitwendige  kracht. 

Indien  er  plotseling  een  kracht  op  dit  systeem  wordt  geplaatst  dan  volgt 
een  verplaatsing; 


X 


f 


-  1 
K  \ 


-  c  t 

-  e  cos(cOj  t) 


7 


De  gedempte  sinusvorm  heeft  een  cirkel frequentie : 


CO 


1 


K 


^  M 


^  C 
2M 


V 

7 


Zonder  demping  blijft  de  amplitude  gelijk  en  is  de  cirkelfrequentie : 


(Oq  = 


N 


K 

~M 


Voor  een  kritische  waarde  van  de  demping  C^rit  degenereert  de  oplossing  tot 
een  uitdempende  e-macht  (  (Oi=0  )  : 
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Het  is  gebruikelijk  om  de  damping  uit  te  drukken  als  een  fractie  P  van  deze 
kritische  damping; 


P  = 


c 

2Ma>o 


C(Oo 

2K 


Hiarmee  wordt  da  cirkalfraquantia  ©i: 


COj 


Waaruit  volgt  dat  ©i  niet  vaal  varandert. 

Da  waarda  van  P  kan  ook  warden  afgeleid  uit  de  verhouding  van  de  amplitudes 
van  opvolgende  maxima  (of  minima)  x^.j^  an  x^  in  de  verplaatsing  (snelheid  of 
versnelling) : 


P  = 


2n 


1  X. 

In(-L) 


^,-1 


Daze  formula  is  gebruikt  om  de  damping  te  berekenen  uit  de  berekende 
verplaatsingen  of  versnell ingen. 

Voor  een  constructie  gelden  vergelijkbare  formulas  maar  de  trilvormen  zijn 
dan  veal  complexer  an  ook  zijn  ar  maar  aiganfraquantias  waarin  da 
constructia  kan  trillan.  Indian  de  eindige  elementen  methode  wordt  gebruikt 
om  een  constructie  door  te  rekenen  dan  zijn  M,  C,  an  K  matrices  an  x 
bestaat  uit  aen  vactor  van  varplaatsingan  van  punten  van  de  constructie  en 
F  bestaat  uit  een  vector  van  uitwendige  krachten. 

Het  is  in  dit  geval  niet  eenvoudig  om  eigenfrequenties,  trilvormen  en 
damping  te  berekenen.  Uitgaande  van  een  sinusvormige  verplaatsing  van  een 
dekveld  met  voornoemde  dempingsveren  met  constante  K  en  knoopuntsmassa  Mq 
volgt  een  benadering  voor  P ; 


P  = 


n  K 

TmJJp 


waarin  fg  de  eigenfrequentie  van  een  plaatveld  en  Mq  de  knooppuntsmassa  {in 
ons  geval  85  kg)  en  N  het  aantal  elementen  in  het  plaatveld. 

Voor  de  zijvelden  geldt  N=8,  fo=33  Hz  en  de  benaderde  damping  is  .028. 

Voor  het  middenveld  geldt  N=12,  fo=14  Hz  en  demping  .03. 

Deze  waarden  komen  redelijk  overeen  met  de  dempingswaarden  van  4  resp.  3  % 
die  uit  de  computerberekeningen  kwamen.  De  aangenomen  sinusvormige 
verplaatsing  past  overigens  het  beste  bij  de  trilvorm  met  de  frequentie  van 
14  Hz. 

Dezelfde  dempingsveren  zijn  gebruikt  in  som  15. 
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an  actual  charge  geometry.  Cable  wound  shock  mounts  were  installed.  Shock  response  results  were  compared  with  those  as  obtained  earlier  for  the 
same  geometry.  Stiffness  of  the  middle  deck  was  changed  in  order  to  study  the  effect  of  nearly  coinciding  eigen  hrequencies  and  the  influence  of 
damping  was  studied. 
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